arXiv:1507.08550vl [astro-ph.GA] 30Jul2015 


Modelling dust extinction in the Magellanic Clouds 


Alberto Zonca, Silvia Casu, Giacomo Mulas, Giambattista Aresu 

INAF-Osservatorio Astronomico di Cagliari, Via della Scienza 5, 1-09047 Selargius, Italy 

azoncaOoa-cagliari.inaf.it,silviaSoa-cagliari.inaf.it,gmulas0oa-cagliari.inaf.it, 

garesuSoa-cagliari.inaf.it 

Gesare Gecchi-Pestellini 

INAF-Osservatorio Astronomico di Palermo, P.za Parlamento 1, 1-90134 Palermo, Italy 

cecchi-pestelliniQastropa.inaf.it 
Accepted. Received; in original form 


ABSTRACT 

We model the extinction profiles observed in the Small and Large Magellanic clouds with 
a synthetic population of dust grains consisting by core-mantle particles and a collection of 
free-flying polycyclic aromatic hydrocarbons. All different flavors of the extinction curves ob¬ 
served in the Magellanic Clouds can be described by the present model, that has been previously 
(successfully) applied to a large sample of diffuse and translucent lines of sight in the Milky Way. 
We find that in the Magellanic Clouds the extinction produced by classic grains is generally 
larger than absorption by polycyclic aromatic hydrocarbons. Within this model, the non-linear 
far-UV rise is accounted for by polycyclic aromatic hydrocarbons, whose presence in turn is al¬ 
ways associated to a gap in the size distribution of classical particles. This hints either a physical 
connection between (e.g., a common cause for) polycyclic aromatic hydrocarbons and the absence 
of middle-sized dust particles, or the need for an additional component in the model, that can 
account for the non-linear far-UV rise without contributing to the UV bump at ~217 nm, e.g., 
nanodiamonds. 


Subject headings: dust, extinction - evolution - galaxies: ISM 

1. Introduction 


Interstellar Extinction Curves (ISECs) are now 
available for many lines of sight in the Milky Way 
Galaxy (MWG), in the Small and Large Mag¬ 
ellanic Clouds (SMC, LMC), in M31, in nearby 
galaxies, and increasingly in objects at high red- 
shift. The ISECs have a great variety of shapes. 
This evidence may seem to suggest that grain 
models must be complex, allowing for different 
carriers for each of the different components of the 
ISECs. However, all the measured extinction pro¬ 
files are clearly members of the same family, and in 
fact, ISECs can be phenomenologically described 
by means of a very limited number of empirical pa¬ 
rameters (Cardelli et al.|19^ Valencic et al. 2004 


Fitzpatrick fc Massa||2007). Such a smooth topo¬ 


logical transition among extinction profiles may be 
naturally interpreted as the response of dust to the 
local physics. 

Dust grains are composed mainly of silicates 
and carbons, as made clear by spectroscopic ev¬ 
idence and depletion data (e.g., Draine 20031. 


They are formed in the envelopes of cool stars 
and in novae and supernovae, and are destroyed 
in shocks, but with different efficiencies, silicates 


being more resistant than carbons (Zhukovska & 


Henning 2013 Bocchio et al. 20141. Even if the 


silicates and carbons can be produced as sepa¬ 
rate populations, the destruction of solid carbons 
to atoms and small molecules will lead to car- 
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bon deposition on silicate cores. According to 
the scenario originally envisaged by [Jones et al.| 
(|1990p, gaseous carbon in diffuse clouds is grad¬ 


ually deposited on silicate cores in hydrogen rich 
sp^ form, and is eventually annealed and darkened 
to graphitic (H-poor) sp^ kind by the interstellar 
radiation field. In other words, dust is not im¬ 
mutable and is expected to evolve in time (e.g., 
Cecchi-Pestellini et al.|[2M0 |. 

Such dust formation and evolution scenario has 
been implemented in a model - the [CM]^ model 
(Cecchi-Pestellini et al. 20081 - successfully ap¬ 


plied to a sample of 329 galactic lines of sights 


(Fitzpatrick & Massa 20071 by Mulas et al. (20131. 


Subsequently, the carbon evolution prescrip¬ 
tions of this scenario were condensed in a set of 
time-dependent equations, describing carbon de¬ 
position onto the silicate cores and its following 
photo-darkening, and matched to the observation- 
ally based inferences of dust physical parameters 


(Cecchi-Pestellini et al. 

2014), closing the circle. 

Most of the ISECs in 

the 

Fitzpatrick & Massa 


(20071 sample correspond to dust mantles that 


restarted growing after being almost completely 
removed, with the leftover fully rehydrogenated, 
by, e.g., occasional shocks. 

A natural extension of this work is its appli¬ 
cation to (and test with) environments markedly 
different from the interstellar medium of the so¬ 
lar neighborhood, in which the life cycle of the 
interstellar medium is known to be very different 
either globally, e.g., much more frequent shocks 
as in starburst galaxies, or locally, because of in¬ 
dependent observational constraints on individ¬ 
ual lines of sight. In this paper we shall apply 
the [CM]^ model to a number of lines of sight in 
the Magellanic Clouds (MCs), whose data have 


been provided by Gordon et al. (2003) and Car- 


tledge et al. (2005). These two irregular dwarf 
galaxies are nearby 50 kpc), gas-rich, actively 
star-forming, and have sub-solar metallicities (e.g., 
Howarth]|20Il |. The SMC average ISEC presents 
a featureless profile with no 217.5 nm bump and a 
steep far ultraviolet rise, while the LMC one shows 
a hint of the bump and a firm far-ultraviolet rise, 
stronger than in the MWG but less sharp than in 
the SMC. 

We describe the specific lines of sight, their 
known physical properties, and how they are mod¬ 
elled in Section the resulting characterisation 


of dust grains along such extinction paths in Sec¬ 
tion 1^ Section contains our discussion and the 
last Section summarises our conclusions. 


2. Extinction observed and modelled 


Gordon et al. 

2003)and 

Gartledge et al. 

(2005 


exploited a collection of newly taken and archive 
data to obtain individual extinction curves for a 
total of 24 lines of sight: 19 for LMC and 5 
for SMC. The list of sight lines reported in Ta¬ 
ble contains also their absolute visual magni¬ 
tude, Ay, the ratio of total to selective extinc¬ 
tion, Rv, and the gas to dust ratio expressed 
as Nh/Eb-v cm“^ mag“^. This quantity has 
been constructed for specific lines of sight using 
data taken from Welty et al. (2012) and references 
therein: the atomic hydrogen column density is 
derived from absorption Lyman-a profiles in all 
cases but Sk-69 279, for which we used data from 
the hydrogen 21 cm emission; the molecular hy¬ 
drogen column density is obtained from FUSE H 2 
absorption measurements when available, or oth¬ 
erwise exploiting the relation linking N-h^/Eb-v 
to Eb-Vj as found by Welty et al. (2012) for the 
MWG. In addition to individual lines of sight, we 
also considered the “average” extinction curves 
given by Gordon et al. (2003) for the SMG bar, the 
LMG 2 Supershell (Ss), and the overall LMG aver¬ 
age. These average extinction curves were treated 
exactly as the individual ones, to see if this re¬ 
sulted in a sensible description of representative 
dust properties for the areas they refer to. 


The dust model consists in a distribution of 
core-mantle grains and a mixture of polycyclic 
aromatic hydrocarbons (PAHs). Glassical grains 
are constructed assembling four concentric com¬ 
ponents: a cavity, whose volume is in fixed pro¬ 
portion fy for all core sizes; the silicate shell of 
radius a; the sp^ layer of thickness fsp^w, w be¬ 
ing the mantle thickness; and the sp^ layer of 
thickness fgpsw = (I — fgpi)w. We consider a 
size distribution for the silicate cores, given by a 
power law (a -I- Such distribution allows 

for a gap in particle sizes, so that two popula¬ 
tions of dust grains, big and small ones, may be 
present, each characterised by lower (a_, 6_) and 
upper size (a+, 6+) limits in the distribution. The 
molecular component is represented by a mixture 
of 54 PAHs in the size range 10 — 66 G atoms, 
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Table 1: Relevant properties of the lines of sight toward the Magellanic Clouds in our sample (in italic the 
average ISECs). 


LoS 

Av 

(mag) 

Rv 

(xlO' 

Nn/EB-v 

21 cm “2 mag“i) 

AzVlS 

SMC bar sample 

O .491 2.901 

67.193d 

AzV23 

0.48^ 

2.652 

49 . 293 d 

AzV214 

0.352 

2.402 

41.033d 

AzV398 

0.682 

3.142 

55.633d 

SMC Bar 

— 

2.742 

24 . 55 I' 

AzV456 

SMC wing sample 

0.571 2.191 

10 . 393 d 


LMC 

Average Sample 

Sk-66 19 

0.862 

3.442 

29 . 593 d 

Sk-66 88 

1.032 

3.672 

20.713d 

Sk-67 2 

0.561 

3.751 

18.553d 

Sk-68 23 

1.042 

3.352 

6.033d 

Sk-68 26 

O .641 

3.451 

20 . 773 d 

Sk-68 129 

O .571 

3.371 

32 . 733 d 

Sk-69 108 

0.982 

3.152 

10.883d 

Sk-69 206 

0.962 

3.682 

35 . 333 d 

Sk-69 210 

1.362 

3.322 

25.753d 

Sk-69 213 

0.632 

3.962 

22 . 553 d 

LMC 

— 

3.412 

12 .881’ 

Sk-68 140 

LMC2 Ss sample 

0.671 3.341 

28.773d 

Sk-68 155 

0.561 

2 . 8 II 

21.823d 

Sk-69 228 

O .531 

3.541 

28.503d 

Sk-69 256 

0.112 

0.642 

15.673d 

Sk-69 265 

0.322 

1.682 

24 . 433 d 

Sk-69 270 

0.442 

2.342 

18.993d 

Sk-69 279 

O .751 

3.541 

21.80iid 

Sk-69 280 

0.562 

3.122 

23.033d 

Sk-70 116 

0.652 

3.412 

16.823d 

LMC2 Ss 

— 

2.762 

21 .881’ 


Cartledge et al. (2005); ^ Gordon et al. (2003); atomic hydrogen from Lyman a; molecular hydrogen 


from FUSE] molecular hydrogen estimate (see text); ® atomic hydrogen from 21 cm; ^ Welty et al. (2012) 
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in the charge states 0, ±1, and +2 (Mailed et 


al. 2007). The equations describing the model 


have been reported in several papers (e.g., Cecchi- 


Pestellini et al. 2008 Zonca et al. 2011), and are 


not repeated here. The model produces a theoret¬ 
ical extinction curve as a function of its parame¬ 
ters. A least-squares procedure (Levenberg [1944 


Markwardt|2009 ) yields the best-fitting model pa¬ 


rameter values for each line of sight. The vari¬ 
ance-covariance matrix for the best-fitting param¬ 
eters, measuring how tightly they are constrained 
by the fit, is obtained by synthetic statistics, by 
repeating the fit multiple times on data perturbed 
according to their errors. The details of this proce¬ 
dure, and the modifications in the MPFIdj^ Lev- 
enberg-Marquardt implementation to cope with 
underdetermination of some parameters, are de¬ 
scribed in Mulas et al. (2013). 

The evolutionary model of carbon mantles is 
detailed in Cecchi-Pestellini et al. (2014). The 
carbon cycle in and out of dust is described by a 
set of ordinary differential equations that follow 
the evolution in time of the available gas-phase 
carbon, the deposited polymeric sp^ carbon , 
and its radiation-annealed sp^ counterpart. The 
main free parameters are the photo-darkening 
time — the time scale for graphitization (e.g.. 


Cecchi-Pestellini et al. 2010) — and the initial 


conditions. The differential equations thus consti¬ 
tute a one-parameter system: solutions starting 
from the same initial conditions, but with differ¬ 
ent photo-darkening times, have no other points 
in common but the starting point and the asymp¬ 
totic limit. Of course, the evolutionary model also 
has a parametric dependence on the assumed en¬ 
vironmental conditions, i.e. the elemental abun¬ 
dances and so on. However, this dependence is 
rather mild, and the solutions of the model re¬ 
tain their overall behaviour with minor differences 


for any sensible environmental conditions (Cecchi- 
Pestellini et al.pOld ). 


3. Results 

Fitted ISECs are shown in Figs. (SMC) and 
(LMC), while the inferred dust parameters are 
reported in Tables [^(classical dust) an d||(PAHs). 
In the same tables we also show the results for the 
average extinction profiles. 


^http://purl.com/net/mpfit 



Fig. 1.— Fits to the normalized extinction curves 
for five lines of sight in the Small Magellanic 
Cloud. Solid black lines: global fit. Dashed lines: 
classical dust contribution; in blue the contribu¬ 
tion of “large” dust grains, in red “small” dust 
grains, in black their sum; only the black curve 
is shown when the gap is very small or nonexis¬ 
tent. Dot-dashed lines: PAHs. For lines of sight 
in the SMC bar, we also show, in solid purple, the 


average” extinction curve for the SMC bar (Gor- 
don et al.||2003). The lightly shaded areas are the 


observational error ranges (see text). 
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Fig. 2.— Fits to the normalized extinction curves for lines of sight in the Large Magellanic Cloud. Solid 
black lines: global fit. Dashed lines: classical dust contribution; in blue the contribution of “large” dust 
grains, in red “small” dust grains, in black their sum; only the black curve is shown when the gap is very 
small or nonexistent. Dot-dashed lines: PAHs. In solid purple, the “average” extinction curve for the LMC 
( Gordon et al.|[2003 ).; dot-dashed lines: PAHs. The lightly shaded areas are the observational error ranges 
(see text). 
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Fig. 2.— Fits to the normalized extinction curves for lines of sight in the Large Magellanic Cloud LMC2 
Ss. Solid black lines: global fit. Dashed lines: classical dust contribution; in blue the contribution of “large” 
dust grains, in red “small” dust grains, in black their sum; only the black curve is shown when the gap is 
very small or nonexistent Dot-dashed lines: PAHs. In solid purple, the “average” extinction curve for the 
LMC (Gordon et al. 2003).; dot-dashed lines: PAHs. The lightly shaded areas are the observational error 
ranges (see text). 
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Table 2: Classical dust component parameters and Carbon column density locked in PAHs for each line of sight. 
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We adopt si = 0.000001, £2 = 0.01, £3 = 0.0001 
Fit solution for AzV23 without PAHs 
Fit solution for AzV23 including PAHs 

Fit solution for Sk—69 256 with similar contributions by classical dust and PAHs 
Fit solution for Sk-69 256 without PAHs 



Table 3: Mean properties of the PAH mixture (in italic the average ISECs). 


LoS 

N^^^/Av 

(0)/[C] 

^q/[C] 


10^^ (cm“^ mag“^) 

(e-) 

(e-) 


SMC bar sample 


AzVlS 

2.8 (0.4) 

-0.053 (0.012) 

0.051 (0.013) 

AzV214 

1.9 (0.4) 

-0.007 (0.021) 

0.073 (0.010) 

AzV23(“) 

3.72 (0.21) 

0.007 (0.012) 

0.075 (0.004) 

AzV398 

2.6 (0.3) 

-0.010 (0.017) 

0.068 (0.004) 

SMC Bar 

3.4 (0.3) 

0.006 (0.022) 

0.072 (0.010) 


SMC wing sample 


AzV456 

4.40 (0.19) 

0.060 (0.008) 

0.045 (0.008) 


LMC average sample 


Sk-66 19 

1.13 (0.17) 

0.015 (0.020) 

0.063 (0.011) 

Sk-66 88 

1.06 (0.19) 

-0.003 (0.013) 

0.064 (0.006) 

Sk-67 2 

3.43 (0.04) 

0.025 (0.001) 

0.051 (0.001) 

Sk-68 129 

1.6 (0.3) 

-0.012 (0.011) 

0.071 (0.005) 

Sk-68 23 

1.40 (0.08) 

-0.001 (0.009) 

0.063 (0.005) 

Sk-68 26 

1.38 (0.10) 

0.000 (0.010) 

0.056 (0.006) 

Sk-69 108 

1.72 (0.22) 

0.004 (0.016) 

0.055 (0.009) 

Sk-69 206 

0.75 (0.09) 

0.015 (0.008) 

0.056 (0.007) 

Sk-69 210 

1.36 (0.09) 

0.026 (0.006) 

0.050 (0.006) 

Sk-69 213 

2.43 (0.26) 

0.034 (0.010) 

0.065 (0.006) 

LMC 

2.6 (0.6) 

0.019 (0.022) 

0.068 (0.010) 


LMC2 Ss sample 


Sk-68 140 

2.18 (0.17) 

-0.030 (0.010) 

0.066 (0.007) 

Sk-68 155 

1.43 (0.16) 

-0.009 (0.013) 

0.076 (0.005) 

Sk-69 228 

0.53 (0.06) 

-0.013 (0.021) 

0.069 (0.012) 

Sk-69 256('’) 

9.5 (1.8) 

0.05 (0.06) 

0.07 (0.03) 

Sk-69 265 

1.8 (0.4) 

-0.005 (0.007) 

0.061 (0.007) 

Sk-69 270 

4.06 (0.26) 

0.007 (0.009) 

0.073 (0.007) 

Sk-69 279 

0.64 (0.17) 

0.010 (0.018) 

0.065 (0.006) 

Sk-69 280 

2.24 (0.23) 

-0.003 (0.010) 

0.055 (0.007) 

Sk-70 116 

0.94 (0.10) 

-0.006 (0.016) 

0.074 (0.007) 

LMC2 Ss 

2.7 (0.6) 

0.008 (0.020) 

0.068 (0.010) 


We adopt ei = 0.001, £2 = 0.1 

Fit solution for AzV23 including PAHs 

Fit solution for Sk-69 256 with similar contributions by classical dust and PAHs 



Table 4: Total Si and C abundances for the lines of sight plotted in Figs. and All the abundances are 
espressed in ppM. 


LoS 

Si/H (ppM) 


C/H (ppM) 


Si/C 


mantle 

PAH 

total 




SMC bar sample 



AzVlS 

4.8 (1.0) 

2.7 (1.3) 

12.1 (1.6) 

14.8 (1.6) 

0.3 (0.1) 

AzV214 

7.6 (1.2) 

2.0 (0.4) 

11.3 (2.2) 

13.3 (2.5) 

0.6 (0.2) 

AzV23 

6.8 (< e) 

65.0 (0.3) 

— 

65.0 (0.3) 

0.3 (0.1) 

AzV23 

6.9 (1.4) 

0.7 (1.2) 

20.0 (1.1) 

20.7 (1.2) 

0.1 (< e) 

AzV398 

7.1 (1.7) 

3.5 (1.5) 

14.9 (1.8) 

18.4 (2.4) 

0.4 (0.1) 



SMC wing sample 



AzV456 

17.8 (1.1) 

9 (5) 

93 (4) 

102 (7) 

0.2 (< e) 



LMC Average sample 



Sk-66 19 

13.8 (0.6) 

4.0 (1.5) 

13.1 (1.9) 

17.1 (2.9) 

0.8 (0.1) 

Sk-66 88 

20.5 (0.8) 

10.7 (1.2) 

19 (3) 

29 (4) 

0.7 (0.1) 

Sk-67 2 

18.1 (0.7) 

5.4 (0.2) 

69.4 (0.9) 

74.8 (1.1) 

0.2 (< s) 

Sk-68 129 

12.0 (1.4) 

4.0 (2.3) 

16 (3) 

20 (5) 

0.6 (0.2) 

Sk-68 23 

61.5 (0.8) 

14 (4) 

78 (4) 

91 (8) 

0.7 (0.1) 

Sk-68 26 

18.1 (0.6) 

10 (4) 

22.8 (1.7) 

33 (5) 

0.6 (0.1) 

Sk-69 108 

31.0 (1.1) 

9.2 (2.1) 

50 (6) 

59 (8) 

0.5 (0.1) 

Sk-69 206 

12.0 (0.3) 

5.0 (0.8) 

7.9 (1.0) 

12.9 (1.3) 

0.9 (0.1) 

Sk-69 210 

13.5 (0.5) 

4.7 (0.7) 

17.6 (1.1) 

22.3 (1.0) 

0.6 (< e) 

Sk-69 213 

17.4 (1.9) 

5.3 (1.3) 

43 (5) 

48 (4) 

0.4 (0.1) 



LMC2 Ss sample 



Sk-68 140 

11.5 (1.3) 

7.8 (2.7) 

25.3 (2.0) 

33 (4) 

0.4 (0.1) 

Sk-68 155 

17.2 (2.4) 

3.9 (0.7) 

18.4 (2.0) 

22.2 (1.9) 

0.8 (0.1) 

Sk-69 228 

16.1 (2.2) 

12.7(1.1) 

6.5 (0.7) 

19.3 (1.5) 

0.9 (0.2) 

Sk-69 256('=) 

5.0 (1.0) 

18 (10) 

39 (7) 

57 (13) 

0.1 (< e) 

Sk-69 256(‘') 

10.4 (0.4) 

21(5) 

— 

21 (5) 

0.5 (0.2) 

Sk-69 265 

10.0 (1.4) 

2.0 (0.3) 

12.1 (2.4) 

14.2 (2.4) 

0.7 (0.2) 

Sk-69 270 

13.2 (2.7) 

3.9 (1.1) 

50 (3) 

53.9 (2.9) 

0.3 (0.1) 

Sk-69 279 

23 (4) 

8.7 (1.9) 

10.4 (2.8) 

19 (3) 

1.2 (0.3) 

Sk-69 280 

14.7 (1.1) 

2.0 (1.2) 

30 (3) 

32 (3) 

0.5 (0.1) 

Sk-70 116 

23.3 (1.7) 

5.9 (1.0) 

19.0 (2.1) 

25.0 (2.2) 

0.9 (0.1) 


We adopt e = 0.1 

Fit solution for AzV23 obtained with classical dust only (no PAHs) 

Fit solution for AzV23 including PAHs 

Fit solution for Sk-69 256 with similar contributions by classical dust and PAHs 
Fit solution for Sk-69 256 obtained with classical dust only (no PAHs) 
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As described in Mulas et al. (2013), the fit is 
severely underconstrained, mainly due to degener¬ 
acy in how many different PAH mixtures can add 
up to very nearly the same cumulative cross sec¬ 
tion. As a result, only some of the parameters are 
well determined by the fit, namely the ones defin¬ 
ing the classical dust, and some collective proper¬ 
ties of the PAH mixture. Generally, the respective 
contributions of classic dust and PAHs to the mod¬ 
eled ISECs (as shown in Figs. and are very 
well constrained. This was the case for all ISECs 
we fitted in the MWG and for all but two excep¬ 
tions in the MCs, i.e. AzV23 and Sk-69 256. In 
the latter two exceptions, to add insult to injury, 
there is more than one acceptable solution, with 
qualitatively different shapes of the classic dust 
and PAHs contributions that produce almost the 
same profile. These degenerate fits are shown in 
Fig-i and parameters for both solutions are given 
in Tables §[§ an d[5] This reflects an actual phys¬ 
ical degeneracy in the possible ways to account 
for the observed shapes of these specific extinction 
curves, which could be broken by using additional 
information (e.g., abundance constraints or obser¬ 
vations of aromatic features in emission). AzV23 
has an extremely weak bump, if any, and a corre¬ 
spondingly weak non-linear far-UV rise. This can 
be achieved by a plain monomodal distribution 
of classic dust with negligible PAHs (red curve in 
Fig. 0; alternatively the same ISEG results from 
opening a gap in the grain distribution, thereby 
creating a broad break in the classical dust extinc¬ 
tion which almost perfectly cancels with the bump 
due to PAH absorption. To make evident the effect 
on the extinction curve of a gap in the size distri¬ 
bution of classical dust grains, in the top panel of 
Fig. 1^ we show the ISEG towards Sk-67 2, the best 
fitting extinction curve (which has a gap) and its 
components in black, with superimposed, in blue, 
the extinction curve that would result from ex¬ 
actly the same extinction model but with no gap. 
In the case of Sk-69 256 there is a hint of a bump, 
but due also to the low signal-to-noise ratio a sim¬ 
ilar degeneracy of acceptable solutions occurs, i.e. 
PAHs may or may not be present. The remaining 
ISEGs in the MG sample do not appear to suffer 
from such degeneracy, as their fits unambiguously 
and consistently converge to the same solutions re¬ 
gardless of the starting point of the iterative fitting 
procedure. 



Fig. 3.— Alternative (degenerate) acceptable 
fits for the ISEGs towards AzV23 and Sk-69 256. 
Black solid lines: absolute best fit. Red solid 
lines: alternative fit without PAHs. Dashed black 
lines: classical dust contribution for the best fit. 
Dot-dashed lines: PAHs for the best ht. The 
lightly shaded areas are the observational error 
ranges (see text). 
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Fig. 4.— Effect on the extinction curve of the 
presence (or absence) of a gap in the size dis¬ 
tribution of classical dust grains in the [CM]^ 
model. Top panel: observed (shaded area) and 
fitted (solid black curve) ISEC towards Sk-67 2; 
the dashed line shows the contribution of classi¬ 
cal grains to the extinction. In blue, the total 
extinction (solid) and the contribution of classical 
dust grains (dashed line) resulting from exactly 
the same [CM]^ model with the only difference of 
artificially removing the gap in the size distribu¬ 
tion. Bottom panel: extinction produced by clas¬ 
sical dust grains in the best-fitting [CM]^ model 
(black), in the same model with the size gap ar¬ 
tificially removed (blue), and the contribution of 
dust grains with sizes in the gap (red). 


The well-defined PAH collective properties we 
identify (see Mulas et al. 2013 for details) are 
the total column density of carbon atoms locked 
in PAHs per unit visual magnitude, N^^^/Av 
(or given the gas to dust ratio Nq^^/Nh), and 
the fractions in each charge state summarized as 
an average charge per carbon atom, {Q)/[C\ = 
Y.i ^i^i /Li where fV*, qi, and are col¬ 

umn density, charge, and number of carbon atoms 
for the i—th PAH molecule, respectively. In Ta- 
blej^we also report the charge dispersion, (Tq/[C], 
a measure of the homogeneity of the charging 
within the PAH distribution. The total absolute 
column density of carbon locked in PAHs is not 
an outcome of the fit, since we are dealing with 
normalized extinction, but it can derived through 
the absolute extinction (see Table [^. 

The statistics of the MC sample are relatively 
small (24 lines of sight in total). Nevertheless, 
some trends are already clear. The fitting results 
presented in Table occupy a slightly larger vol¬ 
ume of the parameter space with respect to the 
galactic corresponding values. In particular, while 
for the MWG lines of sight the size distribution of 
classical grains is always bimodal with no excep¬ 
tions (see Table 1 in Mulas et al. 2013), with a gap 
in size Aa = 6_ — a+ between ~ 30 and 200 nm, 
there are a few MC lines of sight for which the size 
distribution shows a negligible gap (i.e. AzV214, 
Sk-69 228 and one of the degenerate solutions for 
AzV23). In general, while there is considerable 
scatter in both, median values of Aa are smaller 
in the MCs with respect to those we previously 
found for the MWG. The same can be seen if one 
compares the fits for the average MWG ISEC to 
the average ISECS for the MCs. 

Another important difference between the two 
samples resides in the predominant presence of 
sp^ carbon mantles in the MCs. The sp^ fraction 
is different from zero only in four cases, namely 
Sk-68 26, Sk-69 256, Sk-69 210, and AzV23 (see 
Table [^, while in the MWG carbonaceous man¬ 
tles tend to be overwhelmingly aromatic. This 
is particularly significant as MCs have generally 
stronger radiation fields than the MWG (up to 
ten times higher, [Israel et M||1993 |. 

Erom Figs. and it is evident that classical 
dust dominates the extinction profile, with some 
exceptions in the bump spectral region. In only 
one ISEC (Sk-69 256) PAHs may be dominant in 
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the whole range, depending on which of the for¬ 
mally acceptable solutions is the real one. 

The PAH charge dispersion decreases with in¬ 
creasing PAH charge (Fig.[^ suggesting that when 
the mixture charge is not negligible all molecular 
species tend to assume the same charge value, as 
characteristic during e.g., photoionization (mod¬ 
ulo the differences in the photo-absorption cross- 
sections). On the other side, large charge disper¬ 
sions are expected when a superposition of differ¬ 
ent environments, with different physical condi¬ 
tions, are crossed by a line of sight. In this case 
the average charges of individual regions merge to 
a median, close to zero, value, but with a large 
spread around it. 


4. Discussion 


Figurej^shows in graphical form the total abun¬ 
dances of C and Si locked in dust and PAHs that 
we obtained in this work for the sightlines to the 
MCs examined here, also listed in Table com¬ 
pared with the corresponding values we previously 
obtained for the MWG (Mulas et al. |2013[ ) and 
with some reference values for the total elemen¬ 
tal abundances. The level of compatibility of the 
amounts of Si and C locked in dust in the MCs 
with abundance constraints is at least as good 
as that for the MWG, with a similar fraction of 
data points exceeding mildly the nominal (aver¬ 
age) budget of available Si, C constraints almost 
always respected, one spectacular failure for LMC, 
one for the SMC and one for the MWG (not visi¬ 
ble, because it falls out of the plot). This can be 
interpreted as a validation of the [CM]^ model for 
the MCs. 

Table reports the column mass densities de¬ 
rived from our fits of the individual extinction 
curves, while Figure shows the correlation of 
the total dust mass column densities versus to¬ 
tal hydrogen mass column density. They can be 
compared with the results of |Roman-Duval et al.| 
(20141, obtained fitting a dust emission model to 


Herschel data. While of course our data points are 
very sparse, compared to all the pixels in the maps 
of the MCs, as used by Roman-Duval et ^ ( 2014| ) , 
they do fall in the same area of the Sdust-S(FItot) 
plane covered by their data, and are therefore fully 
consistent with them. This is rather interesting, 
since the two datasets were obtained using com¬ 


pletely different observations (emission versus ab¬ 
sorption) and dust models. 


In Cecchi-Pestellini et al. (20141 we outlined a 
self-consistent evolutionary scenario for dust in 
the diffuse ISM (DISM) in the MWG. In such 
model, silicate grain particles, being relatively 
long-lived, progressively accrete C atoms in an 
H-rich environment (the DISM), thereby form¬ 
ing aliphatic carbonaceous mantles. The older 
layers of these mantles are progressively photo- 
processed to aromatic, while new aliphatic ma¬ 
terial simultaneously continue to deposit on top 
of them. The deposition of carbon is assumed 
to be proportional to the rate of collisions of C 
atoms onto dust grains, with a sticking factor 
calibrated to match observational constraints on 
CO abundances. Since dust evolution is modelled 
in diffuse and translucent clouds, which are rich 
in atomic H, all freshly deposited carbonaceous 
material is assumed to be highly hydrogenated, 
i. e. aliphatic. Upon UV irradiation, carbona¬ 
ceous material is assumed to be dehydrogenated, 
becoming more and more aromatic. To sim¬ 
plify the model, carbonaceous material, instead of 
having a continuous, smooth transition between 
deeper, older, and thus more aromatic material, 
and freshly deposited, aliphatic material, is rep¬ 
resented as just two layers, the deeper one com¬ 
pletely aromatic, the other completely aliphatic. 
Occasionally, shocks peel away most of this ma¬ 
terial, restarting the cycle and ejecting the frag¬ 
ments of carbonaceous mantles into the gas, where 
they can possibly be a source of PAHs. Shocks 
are not explicitely included in the set of differ¬ 
ential equations giving the time evolution of the 
grain mantles, but are taken into account in terms 
of “instantaneous” episodes, removing most of the 
carbonaceous mantles, rehydrogenating the small 
leftover, and putting back into the gas phase the 
resulting carbon; the time evolution, according to 
the differential equations of the model, is then re¬ 
sumed from these new initial conditions. 


Upon trying to use this kind of evolutionary sce¬ 
nario to interpret the ISECs observed in the MGs, 
as shown in Fig|^ it is apparent that the anneal¬ 
ing time scale must be almost always much longer 
than the average time between mantle-shattering 
events. The annealing scale is defined by the com¬ 
petition between energetic processing by UV ra¬ 
diation (and cosmic rays) and re-hydrogenation 
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Table 5: Mass column densities of silicates, carbonaceous grain mantles in classical dust, and PAHs, in 
Mq/pc^, as inferred from our fits of individual extinction curves. 


LoS 

S (silicates) 

I](sp^ carbon) 

i;(sp^carbon) 

I](PAHs) 

AzVlS 

0.075 (0.023) 

SMC bar sample 
0 (0.0) 

0.0034 (0.0020) 

0.014 (0.004) 

AzV214 

0.063 (0.013) 

0 (0.0) 

0.0013 (0.0004) 

0.0068 (0.0016) 

AzV23“ 

0.083 (0.005) 

0.049 (0.003) 

0.0071 (0.0005) 

0 (0.0) 

AzV23'’ 

0.084 (0.018) 

0 (0.0) 

0.0007 (0.0012) 

0.0180 (0.0015) 

AzV398 

0.12 (0.03) 

0.000003 (0.000016) 

0.0047 (0.0021) 

0.018 (0.004) 

AzV456 

0.066 (0.010) 

SMC wing sample 
0 (0.0) 

0.0028 (0.0015) 

0.025 (0.004) 

Sk-66 19 

0.140 (0.019) 

LMC Average sample 

0.00001 (0.00003) 0.0033 (0.0013) 

0.0097 (0.0019) 

Sk-66 88 

0.163 (0.020) 

0.000001 (0.000004) 

0.0070 (0.0011) 

0.0109 (0.0023) 

Sk-67 2 

0.069 (0.028) 

0 (0.0) 

0.0017 (0.0007) 

0.019 (0.008) 

Sk-68 129 

0.09 (0.04) 

0.00000 (0.0000011) 

0.0025 (0.0017) 

0.009 (0.004) 

Sk-68 23 

0.158 (0.012) 

0.000002 (0.000005) 

0.0028 (0.0008) 

0.0146 (0.0014) 

Sk-68 26 

0.095 (0.024) 

0.0019 (0.0016) 

0.0022 (0.0009) 

0.0088 (0.0023) 

Sk-69 108 

0.144 (0.015) 

0.00002 (0.00010) 

0.0034 (0.0008) 

0.0169 (0.0028) 

Sk-69 206 

0.151 (0.019) 

0.000002 (0.000008) 

0.0052 (0.0010) 

0.0073 (0.0013) 

Sk-69 210 

0.196 (0.017) 

0.0038 (0.0008) 

0.0012 (0.0003) 

0.0187 (0.0019) 

Sk-69 213 

0.086 (0.019) 

0.000005 (0.000022) 

0.0021 (0.0006) 

0.015 (0.003) 

Sk-68 140 

0.09 (0.03) 

LMC2 Ss sample 

0.00000 (0.0000015) 

0.0050 (0.0024) 

0.015 (0.005) 

Sk-68 155 

0.10 (0.04) 

0.0000 (0.0000008) 

0.0019 (0.0007) 

0.0081 (0.0027) 

Sk-69 228 

0.09 (0.04) 

0 (0.0) 

0.0061 (0.0027) 

0.0028 (0.0013) 

Sk-69 256= 

0.018 (0.008) 

0.0024 (0.0028) 

0.0026 (0.0016) 

0.010 (0.004) 

Sk-69 256=* 

0.038 (0.014) 

0.0037 (0.0022) 

0.0019 (0.0009) 

0 (0.0) 

Sk-69 265 

0.064 (0.015) 

0 (0.0) 

0.00107 (0.00024) 

0.0057 (0.0015) 

Sk-69 270 

0.065 (0.018) 

0.00001 (0.00004) 

0.0015 (0.0005) 

0.018 (0.004) 

Sk-69 279 

0.14 (0.05) 

0 (0.0) 

0.0045 (0.0017) 

0.0048 (0.0020) 

Sk-69 280 

0.084 (0.016) 

0 (0.0) 

0.0009 (0.0006) 

0.0126 (0.0026) 

Sk-70 116 

0.102 (0.019) 

0.000001 (0.000003) 

0.0021 (0.0005) 

0.0062 (0.0012) 


Fit solution for AzV23 obtained with classical dust only (no PAHs) 

Fit solution for AzV23 including PAHs 

Fit solution for Sk-69 256 with similar contributions by classical dust and PAHs 
Fit solution for Sk-69 256 obtained with classical dust only (no PAHs) 
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by hot H atomic gas, turning aromatic carbona¬ 
ceous material back into aliphatic. In the MWG, 
re-hydrogenation was deemed to be by and large 
negligible, but this may be not the case in the 
MCs. Carbonaceous mantles in the MCs re¬ 
sult from the fit to be almost always completely 
aliphatic, with only few exceptions. This points to 
significantly higher star-formation activities, and 
associated supernova rates, and/or much longer 
photoprocessing timescales, possibly due to more 
effective re-hydrogenation. This is in agreement 
with the current star formation rates in the MCs 
(~ 0.04 in the SMC, 


grains (Slavin et al. 2004), because larger grains 


Bolatto et al. 2011 and 


0.4 Mq yr ^ in the LMC, Harris & Zaritsky 


20091 that per unit area are an order of magnitude 
higher than in the MWC (~ 1 Mq yr“^ 


Robitaille 


&: Whitney|2010 ). In general, observations of stel¬ 


lar populations in dwarf galaxies suggest that the 
star formation bursts are sporadic, separated by 
millions to billions of years even in isolated sys¬ 
tems (e.g., Weisz et al.||2008 |. According to their 
star-formation histories (Harris & Zaritsky 2004 


2009[ Rezaei et al. 2014), the MCs appear to be 
typical members of this class of galaxies. MC-type 
ISECs may be then prototypical of dwarf galaxies 
of low metallicity, with high gas content (see e.g., 
Cecchi-Pestellini et al.|[^14b ). 

In the MWC DISM the photo-darkening time 
is T^d ~ 10^/x yr ( lida et al.|l984 Mennella et al. 


2001 ), X being local fluctuations around the aver¬ 


age intensity of the interstellar UV radiation field. 
Such time should be compared with lifetime of car¬ 
bon mantles against sputtering (whose estimate is 


v-l 


frequency of~2xl0 “'yr 
derive that dust mantles in the MWC DISM have 
enough time to be converted in aromatic structure. 


Cecchi-Pestellini et al. (2014) show that the evo¬ 


lutionary age of ISECs peaks at times of few Myr 
(assuming typical diffuse interstellar conditions) 
comparable to the probable average age of diffuse 
clouds. Translating these considerations to the 
DISM in the MCs, we may conclude that the car¬ 
bon recycling time must be much shorter than the 
local photo-darkening time, Tpd < 10^ yr (y = 10, 
Israel et al. 1993). However, such result should 


be taken with care, as many effects are concurring 
in dust destruction and formation processes. For 
instance, small grains are more susceptible to de¬ 
struction by sputtering in shock waves than large 


decouple from the gas and are therefore less ex¬ 
posed to ions. In strong shocks, grain-grain inter¬ 
action may lead to shattering and thus to the cre¬ 
ation of smaller grains ( Asano et al.||2013 ), which 
are then more likely to be sputtered away. Hence, 
the timescale of dust destruction may be inversely 
proportional to the abundance of dust, since the 
rate of collisions is proportional to the grain num¬ 
ber density (e.g., Mattsson et al. |Dwek et al. 


2014). 


Although PAHs are clearly less abundant in the 
MCs than in the MWC (see Table |^, a potential 
problem in our description of the carbon cycle in 
galaxies is posed by the the sharp decline of the 
source of aromatic material in the DISM, since 
most of dust mantles are destroyed before they 
experience a substantial annealing. PAH forma¬ 
tion rates in the cold, carbon-rich winds of evolved 
stars are too slow, the injection time being ap¬ 
proximately 2 Cy. Recent Spitzer surveys in SMC 
also revised down the contribution of ACB stars to 
dust production (e.g., Boyer et al. 2012). With¬ 
out a net positive contribution from supernovae 
to the dust budget, this suggests that dust must 
grow in the DISM or be formed by another yet un¬ 


known mechanism (see Matsuura et al. 2011[ and 


McKee 2011). Recent Herschel photometric and 


spectroscopic observations of the supernova 1987A 


by Matsuura et al. (2015) may suggest that super¬ 


novae can be an important source of dust in the 
interstellar medium, from local to high-redshift 
galaxies. However, the point is controversial as it 
depends on the efficiency of supernova remnants in 


Draine 

1987 

), we destroyng interstellar grains. 

Temim et al. 

(2015 


determined dust lifetimes in the MCs exploiting 
the available data of supernova remnants for these 
galaxies. These authors inferred dust lifetimes in 
the MCs significantly shorter then in the MWC, 
as the cumulative dust production rates by ACB 
stars and supernovae are one order of magnitude 
lower than dust destruction in interstellar shocks 
produced by the expansion of supernova remnants 
Conversely, we are left with the puzzle of what 
becomes of the aliphatic fragments released into 
the DISM of the MCs by the shattering events. 
Duley & Williams (1984) proposed that they 
should be rapidly evaporated into polyyines, and 
thereafter quickly photodissociated under DISM 
conditions. They would therefore have a neg- 
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ligible effect on the extinction. If instead a 
significant population of aliphatic carbonaceous 
nanoparticles were able to survive for a signifi¬ 
cant time, thereby achieving non-negligible abun¬ 
dances, they should be visible both in extinction, 
exhibiting a* a features contributing mainly to 
the non-linear far-UV rise, and given their small 
sizes possibly producing aliphatic C-H emission. 
Indeed, some emission in the aliphatic C-H stretch 


at 3.4 fim is detected toward the MCs (Boulanger 
et al.|[2MT Mori et ^|2012 1. 


Classical dust extinction is generally larger than 
PAH extinction with a few exceptions. In these 
latter cases, the PAH contribution dominates 
only in the bump region. In all these relatively 
“PAH-rich” cases, the classical size distribution 
is strongly bimodal, with a gap ranging from 40 
to 130 nm, with a tendency to increase the non 
linearity of the far-UV rise with increasing size 
breaks. All the remaining ISECs are reproduced 
with classic grain size distributions that tend to 
have markedly lower gaps (see Table [^. 

In Mulas et al.| (2013) we concluded that, in 
the framework of [CM]^ model, bumpless extinc¬ 
tion curves are not necessarily lines of sight devoid 
of PAHs, but rather that the tt* ^ tt PAH ab¬ 
sorption can be very effectively masked by a large 
gap between the maximum size limit of small dust 
grains and the minimum size limit of large ones. 
MWG data were not conclusive since MC-type 
ISECs are rare in our Galaxy. The present results 
suggest a different interpretation based on the fact 
that some ISECs show a more or less weak non¬ 
linear far-UV rise without showing an appreciable 
bump. This is particularly evident in the SMC 
sub sample. In the LMC the most evident trends 
are small bumps with linear far-UV rises. Within 
the current model dust populations, for the SMC- 
type ISEC the only possibility is the introduction 
of PAHs that, through a* <— a absorptions, may 
produce such rise. However, as a by product, PAH 
molecules exhibit tt* •<— tt transitions who gener¬ 
ate the bump. Thus, the only viable solution, for 
bumpless ISECs, is the removal of mid-size grains, 
whose extinction power falls in the bump region, 
from the grain size distribution. 

This may be interpreted in two alternative 
ways. The first one is that we must be missing 
an dominated carbon component, whose sizes 
are intermediate between small solids and macro¬ 


molecules. This is also hinted by the numerical 
tendency during the fitting procedure to go below 
the lowest limit of 5 nm in particle size, that is 
systematic for all the lines of sight in the MWG 
and in the MGs. Such small particles might be 
carbon mantle fragments, possibly produced dur¬ 
ing some destructive events (as discussed in the 
previous paragraph), such as e. g. nanodiamonds 


( Rai fc Rastogi|[MTo| 2012 Rai fc Botet||2014 |. 


The second interpretation is instead that in 
our model we are missing no important compo¬ 
nent contributing to extinction, and hence the 
fact that PAHs are abundant only if middle-sized 
dust grains are absent is real, and must have 
a sound physical explanation. One might ar¬ 
gue that if an interstellar cloud, and the dust 
therein, underwent some shocks in its history, they 
may have shattered middle-sized grains, produc¬ 
ing fragments (small grains) and releasing PAHs 
from the disrupted carbonaceous mantles. This 
would be a consistent explanation for the MWG, 
where carbonaceous mantles are almost entirely 
aromatic, so that their disruption would indeed 
produce PAH-like fragments. But in the MCs 
carbonaceous mantles result to be almost entirely 
aliphatic, so their destruction can hardly be a vi¬ 
able source of PAHs. 


In the case of LMC-type ISECs the results 
of the fitting procedure are more straightforward 
than for the SMC-type, bumpless ones: in this 
galaxy PAHs are present but are a minor com¬ 
ponent, their presence being only evident in the 
bump spectral region. The possible extra aliphatic 
’’mesoscopic” component is marginal in this case. 
Conversely, the problem of the lack of a suit¬ 
able source of aromatic material is therefore more 
pressing for LMC than SMC. 


5. Conclusions 


We here concisely summarise the main conclu¬ 
sions of this work. 


The [CM]^ dust model can very precisely 
match all the ISECs observed in the MCs, 
yielding estimates for the total dust mass 
column density versus total hydrogen col¬ 
umn density that are fully consistent with 
those obtained by |Roman-Duval et al.| 
(2014) in a completely different, indepen¬ 


dent way. 
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Abundance of elements locked up in [CM]^ 
models fitting the ISECs in the MCs are 
at least as compatible with abundance con¬ 


straints as those for the MWG (Mulas et al. 


2013). 


• When ISECs have vanishingly weak bump 
and far-UV non-linear rise, they can be 
equally well fitted either by classical dust 
alone or by classical dust and PAHs. In such 
cases models without PAHs use much more 
Si and C atoms, which may put them at odds 
with abundance constraints. 


Carbonaceous mantles in the MCs resulting 
from the [CM]^ model fits are overwhelm¬ 
ingly aliphatic, whereas those in the MWC 
are predominantly aromatic. Together with 
the observational constraint that radiation 
fields in MCs are typically 10 times stronger 
than the average in the MWC (Israel et 
1993), in the [CM]^ model this implies 


al. 


that mantle-shattering events must be 3> 10 
times more frequent in the MCs than in the 
MWC, and/or much more efficient. Desorp¬ 
tion from such mantles cannot be a signifi¬ 
cant source of PAHs in the MCs. 
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Fig. 5.— The relation between the mean charge 
of the PAH collection and its dispersion for each 
line of sight. 


Fig. 6.— Total abundances of Si and C locked 
in dust and PAHs in ppM. The circle represents 
the average MWG ISEC. The boxes refer to obser¬ 
vational constraints, values consistent with them 
being the ones inside each box: solar taken in |As-| 


plund et a. (2009 

I, Mathis (2000|, B, F, and G 

MWG stars from 

Zubko et al. 

(20041, and SMC 

and LMG from Howarth (2011 

. Gray points are 

derived for ISECs in the MWG ( 

Mulas et al. 20131, 

dark blue points are ISECs in the EMC “average” 
sample, light blue points in the LMC2 Ss, and red 
points in the SMC. 
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Fig. 7.— Correlation plot of the mass column 
densities of classical dust (i.e. excluding PAHs) 
as inferred from our fits of individual extinction 
curves, versus the total mass density of H (atomic 
and molecular, in all charge states) along the same 
individual lines of sight. Blue dots are SMC lines 
of sight, red dots LMC ones. 



Fig. 8.— Observationally derived distribution of 
mantle thickness w and normalized sp^ mantle 
carbon fraction, with relative errors for the SMC 
(red asterisks), the LMC (blue asterisks), and the 
MWG (pale grey symbols). The black symbol 
marked “MW” represents the values correspond¬ 
ing to the average galactic ISEC. Values corre¬ 
sponding to average ISECs in the MCs are also 
indicated. MWG data are taken from IMulas et al.1 

(pMl. 
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